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ABSTRACT 

Aims. This paper studies the connection between the relativistic number density of galaxies down the past light cone in a Friedmann- 
Lemaitre-Robertson-Walker spacetime with non-vanishing cosmological constant and the galaxy luminosity function (LF) data. It 
extends the redshift range of previous results presented in Albani et al., where the galaxy distribution was studied out to z = 1. 
Observational inhomogeneities were detected at this range. This research also searches for LF evolution in the context of the frame- 
work advanced by Ribeiro and Stoeger, further developing the theory linking relativistic cosmology theory and LF data. 
Methods. Selection functions are obtained using the Schechter parameters and redshift parametrization of the galaxy luminosity func- 
tions obtained from an I-band selected dataset of the FORS Deep Field galaxy survey in the redshift range 0.5 < z < 5.0 for its blue 
bands and 0.75 < z < 3.0 for its red ones. Differential number counts, densities and other related observables are obtained, and then 
used with the calculated selection functions to study the empirical radial distribution of the galaxies in a fully relativistic framework. 
Results. The redshift range of the dataset used in this work, which is up to five times larger than the one used in previous studies, 
shows an increased relevance of the relativistic effects of expansion when compared to the evolution of the LF at the higher redshifts. 
The results also agree with the preliminary ones presented in Albani et al., suggesting a power-law behavior of relativistic densities at 
high redshifts when they are defined in terms of the luminosity distance. 

Key words. Cosmology: theory - Galaxies: distances and redshifts ~ Galaxies: luminosity function - large-scale structure of Universe 



1. Introduction 

The galaxy volume number density, that is, the number of galax- 
ies enclosed in a given volume, is a very important quantity in 
cosmology. It gives information about the density of the mass- 
energy in the universe and its evolution, allowing us to test var- 
ious observational features of cosmological models like, for in- 
stance, the observational inhomogeneities of the number counts 
data down our past light cone (Albani et al. 2007). This im- 
portant quantity can be observationally determined by a careful 
analysis of data from galaxy redshift surveys, as has been done in 
a systematic way by several authors who derived the galaxy lu- 
minosity function (LF) <p from such surveys. The LF of galaxies 
is a number density per unit of luminosity (Peacock 1999). The 
LF stemming from the galaxy distribution in a given dataseQ is 
commonly fitted using the profile due to Schechter (1976), 



(f)(1) = f l" e- 



(1) 



where I = L/L., L being the observed luminosity, L, the lumi- 
nosity scale parameter, (p* a normalization parameter and a the 
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' In this paper the term LF refers to the observationally determined 
quantity. It must be noted, however, that the LF does not depend only 
on observations, but also on an assumed cosmology. Thus, although it 
is a way of presenting observations, it also contains theoretical infor- 
mation. It could be considered to be fully observational ;/ the assumed 
cosmology is observationally well-substantiated. 



faint-end slope parameter. Various papers on the LF have used 
a variety of datasets in different wavelength regions: Lin et al. 
(1999), Fried et al. (2001), Blanton et al. (2003), Pozzetti et 
al. (2003), Bell et al. (2003), Norman et al. (2004), Wilmer et 
al. (2006), Ly et al. (2007), and Tzanavaris & Georgantopoulos 
(2008) obtained Schechter-type LF parameters for galaxies with 
redshift values out to z ~ 1.5. Poll et al. (2003), and Rudnick et 
al. (2003) did the same for galaxies with redshift values out to 
z ~ 3. More recently Bouwens et al. (2007), and Gabasch et al. 
(2008) obtained those parameters for galaxies with redshift val- 
ues out to z ~ 6. It was found that the LF evolves with redshift 
in all galaxy redshift surveys. 

Using the LF data to study observational features of cosmo- 
logical models requires, however, some model linking relativis- 
tic cosmology number count theory with LF astronomical data 
and practice. One approach to this theoretical link was developed 
by Ribeiro and Stoeger (2003; hereafter RS03), who started from 
very general relativistic considerations and then specialized their 
theoretical results, applying them to the LF data provided by the 
CNOC2 redshift survey (Lin et al. 1999). This provides an exam- 
ple of how such a link is established, as well as facilitating vari- 
ous consistency tests between the data and the expected number 
count predictions from the assumed cosmological model. In a se- 
quel Albani et al. (2007; hereafter A07) used the same CNOC2 
dataset to actually obtain and analyze various types of observa- 
tional relativistic densities, that is, the ones defined on the ob- 
server's past light cone0 



- Usually a density is just a density, and the adjective "relativistic" 
does not apply. However, in this paper this expression has to do with 
the different ways the density can be defined, depending on the differ- 
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The aim of this work is to further study the observed num- 
ber densities of galaxies down the observer's past light cone. We 
want to focus on observational inhomogeneities of the number 
count data, analyzing these with differential density and inte- 
gral density measures in order to separate out the effects of cos- 
mic expansion from the effect of LF evolution. As has been dis- 
cussed in detail in A07 and in Rangel Lemos & Ribeiro (2008), 
these measures are capable of revealing observational inhomo- 
geneities even in spatially homogeneous cosmological models. 
This is because spatial homogeneity is defined on constant time 
hypersurfaces, and is therefore a built-in feature of Friedmann 
cosmologies, whereas observational homogeneity is defined by 
measures of mass-energy densities remaining unchanged along 
the past light cone. We also seek to investigate if there is sig- 
nificant evidence of power-law behavior in the observed galaxy 
distribution at high values of redshift, since A07 detected such 
a behavior out to z x 1 . The importance of power-laws comes 
from the fact that they are indicative of self-organized criticality 
in dynamical systems. In other words, scale invariant phenom- 
ena such as power laws are emergent features of a distribution, 
this being a mechanism by which complexity arises in nature 
through simple local interactions. 

The study presented here expands the results presented in 
RS03 and A07, since the model specific equations obtained in 
RS03 were applicable only to LF data whose parameters as- 
sumed the Einstein-de Sitter (EdS) cosmology. And although 
the study presented in AO? was not limited to the EdS model, 
it was limited to the CNOC2 survey, which only probed out to 
z = 0.75 and extrapolated to z = 1 . In addition, although both pa- 
pers made some important theoretical connections between the 
relativistic cosmology equations and the LF data for a specific 
cosmology, they omitted important steps. These can be summa- 
rized as follows. The equations and methodology used in A07 for 
dealing with non-zero cosmological constant models are gener- 
alized, fully presented and discussed in this paper An appendix 
containing a detailed algorithm on how to obtain the key results 
is also included. The role of the completeness function J con- 
necting theory and observations, as advanced in RS03, is made 
explicit here, adding clarity to the results and presenting them 
in a more straightforward manner Finally, a comparison of the 
geometrical versus evolutionary effects, the latter being empir- 
ically obtained from the LF, on the relativistic volume number 
densities is discussed, considering a broad spectral classifica- 
tion of late versus early type galaxies. So, this work further ad- 
vances the treatment of both RS03 and A07 and presents a com- 
prehensive and detailed discussion of how one can proceed in 
applying the very general, model independent, relativistic cos- 
mology equations to specific cosmologies, thus linking any cos- 
mological model to LF data. This is necessary if we want to 
use LF data to observationally test more general cosmological 
models (any more general models will be inhomogeneous and/or 
anisotropic), a task we intend to carry out in forthcoming papers. 
Table [U summarizes the new contributions of this paper as com- 
pared with the results obtained in the previous ones. 

In this paper we assume a Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmology with non-zero cosmological con- 
stant A and derive model independent relativistic equations 
which enable us to solve numerically expressions which include 
LF data. The numerical scheme is presented in detail and ap- 
plied to a galaxy redshift survey dataset that assumes the FLRW 



standard cosmological model in its volume and distance defi- 
nitions: the I-band selected luminosity functions of Gabasch et 
al. (2004 and 2006; hereafter G04 and G06 respectively) for 
the FORS Deep Field (FDF). From the Schechter parameters 
and their redshift parametrization appearing in G04 and G06, 
we calculated the selection functions tfr in eight bandwidths and 
in equally spaced redshift intervals of 0.25 spanning the range 
0.5 < z < 5.0, for the blue bands (G04), and 0.75 < z < 3.0 for 
the red bands (G06). Once the selection functions were calcu- 
lated, we were able to obtain the observational number counts, 
differential number counts and differential and integral densities, 
as defined in RS03 and Ribeiro (2005) and further discussed in 
A07 and Rangel Lemos & Ribeiro (2008). These quantities were 
then used to study observational inhomogeneities in the rela- 
tivistic radial distribution of galaxies belonging to that dataset. 
Therefore, the calculations of this paper go to much higher red- 
shift than the ones presented in A07. The results reached here 
generally agree with theirs, also indicating that observational in- 
homogeneities in the relativistic galaxy distributions can arise 
due to geometrical, null-geodesic effects, even in the context 
of spatially homogeneous universes (Ribeiro 1992, 1995, 2001, 
2005; Rangel Lemos & Ribeiro 2008). 

The plan of the paper is as follows. In ^we briefly present 
the FLRW model, basically to set up notation and write essential 
results. Then we apply the general, model independent, equa- 
tions of relativistic cosmology number-count theory to this uni- 
verse model in order to obtain the differential equations gov- 
erning the evolution of the two basic quantities of our analy- 
sis, namely, the scale factor and the relativistic number counts, 
^describes in detail our numerical approach for solving these 
two differential equations, as well as our procedures for deriv- 
ing all other quantities relevant in our analysis. In doing that we 
also explain how they connect to the the scale factor and the 
number counts. In ^we further develop general equations con- 
necting LF data with relativistic cosmology number-count the- 
ory and obtain the selection functions for the FDF dataset. In ^ 
we study the behavior of radial statistics of the relativistic densi- 
ties obtained in the context of the FLRW model with A and 
out to z = 5.0. We present our conclusions in ^ The appendix 
presents an algorithm detailing how the methodology discussed 
in ^ ^and ^can be applied to a given parametrized LF in 
order to obtain the observational relativistic densities discussed 



2. Standard cosmology with non-zero cosmological 
constant 

2. 1 . The scale factor 

We begin by writing the FLRW line element as follows. 



[d9^ + sin^ 6 d<f>^) 



(2) 



where the time-dependent function S - S(t) is the cosmic scale 
factor, k is the curvature parameter {k - H-1,0, -1) and c is 
the light speed. As is well known, the Einstein's field equa- 
tions with the metric corresponding to this line element yields 
the Friedmann equation, which, if the cosmological constant A 
is included, may be written (e.g. Roos 1994), 



ent measures of distance we use, since in relativistic cosmology a dis- 
tance can be defined by different measures. In this context the adjective 
"cosmological" could also be used. See 33.4l andl 3.5l below. 



^2 ^ ^^Gp„ 



A kc^ 



(3) 
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where p„ is the matter density and we have assumed the usual 
definition for the Hubble parameter, 



H{t) = 



1 dSif) 
Yit) dt 



(4) 



Let us now define the vacuum energy density in terms of the 
cosmological constant. 



Pa = 



A 
8^' 



Since the critical density at the present time is given by. 



PO,c = 



8nG ' 



(5) 



(6) 



where Hq is the Hubble constant, the following relative-to- 
critical density parameter relations hold. 



PO.c 



Pmg PAq 
PO,c PO,c 



(7) 



We have used the zero index to indicate observable quantities at 
the present time. Notice that since A is a constant, then pA = pA„ . 
Considering the definitions (|7]i, we can rewrite the Friedmann 
equation (O at the present time as follows. 



kc^ = //o'5o'(Qo - 1). 



(8) 



In addition, from the law of conservation of energy applied to 
the zero pressure, matter dominated era, we know that. 



Pm ^ S ~ ^ p,n„ BC Sq 

which leads to. 



-3 



Pm 

Pmo 



Si 

53 



(9) 



(10) 



since the matter-density parameter can also be written in terms 
of the critical density as. 



_Pm_ _ SttG 

''HI — ~ 1 ttO P" 

Pc 3//2 



(11) 



We can rewrite equation ([3J as a first order ordinary differential 
equation in terms of the scale factor S it), by using the results in 
equations (l4li- (fTTl i. yielding. 



dS 

-dF-"' 



^maS 



1/2 



The problem we deal with in this paper requires the solution 
of the equation above as well as of another differential equation 
for the cumulative number count (see eq.|20]below). The sim- 
plest approach is to simultaneously solve these two differential 
equations by numerical means. 

2.2. Relativistic number counts 

Let us write the completely general, cosmological, model- 
independent expression derived by Ellis (1971) for the number 
count of cosmological sources, which takes fully into account 
relativistic effects, as follows. 



dN = (dj^dn[n(-k"ua)]p dy. 



(13) 



Here dN is the number of cosmological sources in a volume 
section at a point P down the null cone, n is the number den- 
sity of radiating sources per unit of proper volume in a section 
of a bundle of light rays converging towards the observer and 
subtending a solid angle dQ. at the observer's position, c/^ is the 
area distance of this section from the observer's viewpoint (also 
known as angular diameter distance, observer area distance and 
corrected luminosity distance), m" is the observer's 4-velocity, 
is the tangent vector along the light rays and y is the affine pa- 
rameter distance down the light cone constituting the bundle (see 
RS03, §2. 1, figure 1). The number density, that is, the number of 
cosmological sources per proper volume unit, can be related to 
the matter density p„, by means of. 



Pm 

AL 



(14) 



where, Aig is simply the average galaxy rest mass, dark matter 
included. 

One should point out that the details of the galaxy mass func- 
tion and how it evolves with the redshift are imprinted in the 
LF itself. As a consequence, it will be included in any selection 
functions stemming from the LF. For the empirical purposes of 
this paper this equation is correct to an order of magnitude and 
should be regarded only as an estimate. This means that equa- 
tion (fl4l) enables us to connect the theoretical relativistic quan- 
tities to the LF data and thus include the redshift evolution of 
the mass function empirically (more on this point at the end of 
^4.1l below). However, to actually extract from the LF its implicit 
galaxy mass function requires some assumption about the func- 
tion Aigiz). We shall deal with this problem in a forthcoming 
paper. 

So, for the empirical approach of this paper it is enough to as- 
sume a constant average galaxy rest mass with the working value 
of Aig ~ 1O"M0 based on the estimate by Sparke & Gallagher 
(2000). Hence, if we use equations ( fTOt and ( fTTT i we can rewrite 
equation (fT4l l as. 



SttGM, 



1 



(15) 



The past radial null geodesic in the geometry given by metric 
^ may be written 



dt 
dr 



S 



cVl -kr^ 



(16) 



Since both coordinates t and r are functions of the affine param- 
eter y along the past null cone, we have that 



(12) dt dt dr 



S \ dr 
dy drdy \cVl -kr^l dy' 



(17) 



Assuming now that both source and observer are comoving, then 
u" - c 5{)' and the following results hold. 



• k''Ua = -C k''g()a = -c k° 



dt 



S \ dr 
dy \ Vl -kr^J dy' 



(18) 



In addition, remembering that the area distance is defined by 
means of a relation between the intrinsically measured cross- 
sectional area element dcr of the source and the observed solid 
angle t/^o (ElUs 1971, 2007; Plebanski & Krasinski 2006), we 
have that 



W )2 ^ ^ ^ S'r\d9' + Sin' edcp') ^ 2 
dQ {de^ + sm^Odf) 



(19) 
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We note that for this metric the area distance, a strictly observa- 
tional distance definition, equals the proper distance dp,, a rela- 
tivistic one. Considering equation ([8]l and substituting equations 
(fTSl) . (fTSl l and ( fT9b into equation ( fT3] l. and remembering that here 
c/Qo — An, we are then able to write the number of sources along 
the past light cone in terms of the radial coordinate r and in the 
FLRW model defined in ^2.11 This quantity is given by the fol- 
lowing expression. 



dr ~' 



2GM, 



3. Numerical problem 

The differential equations for the scale factor S and the cumula- 
tive number count A^, equations (fT2l i and ( |20] | respectively, pro- 
vide the basic quantities necessary here, since the expressions 
which actually connect the relativistic theory to LF data are built 
upon them. However, before these two diff'erential equations can 
be solved numerically, some further algebraic manipulations are 
necessary. Next we shall describe these steps and how other 
quantities relevant to our analysis are written in terms of these 
two functions. 

The numerical problem can be summarized as follows. We 
take the radial coordinate r as the independent variable in or- 
der to numerically obtain the functions N(r) and S (r) along the 
past light cone. Then all other quantities of interest are writ- 
ten in terms of these two functions, meaning that numerical re- 
sults for N{r) and S (r) allow us to straightforwardly obtain z (r) 
and dN/dz (r), as well as various cosmological distances, their 
derivatives, observational volumes and number densities. Thus, 
all quantities used in our analysis end up being written in terms 
of the radial coordinate r. 



The two differential equations ( l20b and ( l23T l comprise our 
numerical problem. Solving them simultaneously enables us to 
generate tables for r, S and A^. A computer code using the fourth- 
order Runge-Kutta method is good enough to successfully carry 
out the numerical tasks. The initial conditions ro and A^o are set 
to zero, whereas So can be derived considering that as r — > the 
spacetime is approximately Euclidean, that is, k x Q. This leads, 
from equation ( fTSl l. to ct - -r as well as So - 1- 



(20) 3,2. Differential number count 



The redshift z can be written as 

So 



(24) 



where it is clear that a numerical solution of the scale factor S (r) 
immediately gives us the numerical solution for z{r). We can 
derive the differential number counts dN/dz by means of the ex- 
pression. 



dN_ _ drdN_ 
dz dz dr ' 

with the help of the useful relation 

dr dr dS 
dz dS dz 



(25) 



(26) 



Numerically, since we build all our quantities using N(r) and 
S (r), any derivative with respect to the redshift we wish to eval- 
uate will be similarly written in terms of the derivative of that 
quantity in terms of the radial coordinate. The derivatives in 
equation (|26] | can be taken from definition (l24l i and equation 
(l23T l. enabling us to write. 



3.1. Scale factor 

A simple and straightforward approach to solving the differential 
equations for the scale factor S and the cumulative number count 
A^ is to take the radial coordinate r as the independent variable 
and simultaneously integrate them by numerical means. To do 
so we shall proceed as follows. 

Equation ( fT2l l has its time coordinate implicitly defined in 
the scale factor, therefore we need to rewrite that differential 
equation in such a way that the independent variable becomes 
explicit. Considering equation (jSj, the null geodesic ( fTSl l be- 
comes. 



dr S- 



dz So Ho 



(Qa„)54 - SoH^o - l)S^ + (Q«„5o3)5 



c2-i/o'5o'(no-l)r2 



(27) 



This, together with equation ( l20l l. allows us to write equation 

I as. 



dN /3cQ,„„//o5o' 



dz 



2GM, 



2c2 



r'S 



^(Qa„)54 - SoH^o - 1)52 + (a,„5o')5 



(28) 



dt 
dr 



Ho^S o'^iQ.Q 



ly- 



(21) 3.3. Proper and comoving volumes 



Along the past light cone, we can write the scale factor in terms 
of the radial coordinate as 



(22) 



dS _ dS dt 
dr dt dr 

Thus, we are able to rewrite equation (fl2b in terms of the radial 
coordinate. 



f = -«» 

dr 



(QaJ^^ - 5o'(no - 1)5^ + (Q„,„5o')5 



• H[)^S Q^{Q.o ■ 



l)r2 



(23) 



To find solutions for S(r) we must assume numerical values for 
f2„,(,, Qa„ and Hq. In this paper we adopt the FLRW model with 



Q,„„ = 0.3, Qa„ 



: 0.7 and //q = 70 km s"' Mpc" 



As discussed in RS03, most cosmological densities obtained 
from astronomical observations assume comoving volumes, 
whereas densities derived from theory often assume the local, or 
proper, volumes. The luminosity function, for instance, is nowa- 
days always obtained from galaxy catalogues by assuming the 
comoving volume in its calculation. So, we are only able to com- 
pare those observationally derived parameters with theory if we 
carry out a conversion of volume units. From metric (|2]l it is ob- 
vious that. 



Vl - kr^ 



r^dr sin 9d0d(p = S^dVc 



(29) 



This equation clearly defines the conversion factor between these 
two volume definitions. 
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3.4. Distance measures dido) 

So far we have used only the area distance d,^ as distance defi- 
nition. However, other cosmological distances can, and will, be 
used later. They can be easily obtained from the area distance by 
invoking Etherington's reciprocity law (Etherington 1933; Ellis 
1971,2007), 

d,^{\+zfd,^(\+z)da. (30) 

Here di is the luminosity distance and da is the galaxy-area 
distance. The latter is also known as angular-diameter dis- 
tance, transverse comoving distance or proper-motion distance. 
A fourth distance will also be useful later, the redshift distance 
c/j, defined by the following equation, 

d, = ^. (31) 
no 



Although the reciprocity law is independent of any cosmological 
model, the detailed calculations presented in the previous sec- 
tions are for the FLRW cosmology only. This means that in de- 
riving c/,,, di, da and d, we write the FLRW expression for d,^ and 
derive the others using the general, model independent, expres- 
sion ( l30l l, as well as equation ( |3TI ). Thus, considering equation 
( fT9] ) for d^ and the reciprocity theorem ( l30b , it is straightforward 
to write the other cosmological distances in terms of the scale 
factor, as follows. 





(32) 


Sq r. 


(33) 


Ho\S I 


(34) 



The derivatives of each distance with respect to the redshift 
stiU need to be determined since they will be necessary later. 
Starting with the area distance li^, they can be easily obtained 
from equations (1% and (l24l) . Then it follows that, 



didji) dS dr did^) 



'So 



„, dS 
r + S[ — 

dr 



dz dz dS dr 
Using equation ( l23l l. this expression may be rewritten as. 



(35) 



didn) 
dz 



So 



c^-Ho^SoH^o-iy 



(QaJ54 - 5o'(Qo - 1)^2 + (Q,„„5o')5 



(36) 



The other observational distances can be numerically calcu- 
lated from equation (l35T l if we consider the reciprocity law (l30t . 
Therefore, we have. 



dz 
d(da) 
dz 



= 2(l+z)J, + (l+z)2^, 
dz 



= d,, -H (1 -Hz) 



d(dd 
dz 



(37) 



(38) 



These two equations can also be rewritten in terms of the scale 
factor if we consider equations ( l24l i and (l36l l. yielding. 



dz 



c^-Ho'-SoH^o-iy 



(Qa„)54 -So\^o - 1)52 + {n„„SolS 



(39) 



c^-Ho^SoH^lo-iy 



dz Ho \ (Qa„)54 - SoH^o - 1)52 + (n,„„5o')5 



(40) 



From now on we shall generically use d^ to indicate any ob- 
servational distance, which can be any one of the four cosmo- 
logical distances defined above (/ = a, g, l, z). 

3.5. Differential and integral densities 

The differential density y, gives the rate of growth in number 
counts, or more exactly in their density, as one moves along the 
observational distance c/,. It is defined by (Ribeiro 2005; Albani 
et al. 2007; Rangel Lemos & Ribeiro 2008), 



1 



47r(4)2 d(d,y 



(41) 



whereas the integrated differential density, or simply integral 
density, gives the number of sources per unit of observational 
volume located inside the observer's past light cone out to a dis- 
tance d;. Is is written as. 



where V, is the observational volume. 



V. = -7T(dy 



(42) 



(43) 



These quantities are useful in determining whether or not, 
and within what ranges, a spatially homogeneous cosmological 
model can or cannot be observationally homogeneous as well 
(Ribeiro 1995, 2005; Rangel Lemos & Ribeiro 2008). This is 
because these densities behave very dififerently depending on 
the distance measure used in their definitions - that is, they 
show a strong dependence on the cosmological distance adopted. 
Therefore, as discussed in Rangel Lemos & Ribeiro (2008), 
these measures are the ones we employ in this paper, because 
they are capable of probing the possible observational inhomo- 
geneity of the number counts. 

From a numerical viewpoint it is preferable to write the den- 
sities given by equations iAli and (l42l l in terms of the redshift. 
Thus, the diff'erential density (l4Tl l may be written as. 



dN 



d(di) 



^'^^r''^^'^ dz 



(44) 



A final point still needs to be discussed. The integral density 
(l42l i is a result of integrating % over an observational volume. 
The simplest way of numerically deriving it is shown in what 
follows. Let us differentiate y* in terms of the volume, so that. 



dV, 



(45) 



Considering equation (l44l) . this result leads to the following ex- 
pression. 



diy*V,) _ diy*V)dVi _ dVj _ dN 
dz dVi dz dz dz 

Similarly, it is simple to show that, 
d{y*V,) dN 



dr 



dr 



(46) 



(47) 
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Finally, from these two equations above, as well as from the def- 
initions of y, and y*, it is easy to conclude that the following 
expression holds. 



N 



(48) 



So, the numerical solution of equation ( l20l l together with the 
numerical determination of all distances, as given by equations 
(fT9] l. (|32] |. ( l33T l. (|34] |. allow us to calculate the volume (|43]l and 
evaluate y*. These results fully determine the numerical problem 
for the cosmological model under study. 

In conclusion, once N{r) and S{r) are calculated along the 
past light cone in terms of the radial coordinate r, all other quan- 
tities are straightforwardly obtained with the same functional de- 
pendence: z (r), dN/dz (r), d, (r), d(di)/dz (r), y,- (r), V,- (r) and 
7* (r)- 



4. Observational quantities 

We shall now specialize the equations above to obtain their ob- 
servational counterparts based on the LF data from a specific 
galaxy catalog. As a direct consequence, we show how these re- 
sults can be used to test the consistency of the number count 
theory detailed above. 



4.1. General equations 

Generally speaking, we shall assume that an observational quan- 
tity [T]„i„ can be related to its theoretical counterpart T by means 
of a completeness function J, such that. 



[TU,. = JT. 



(49) 



RS03 showed that such a completeness function can be obtained 
by relating the selection function tp, which gives the number of 
galaxies with luminosity above a given threshold in a given co- 
moving volume, to the radial number density ric (z) in comov- 
ing volume units - that is, the number of sources in a given co- 
moving volume. As akeady mentioned, it is common practice 
to adopt this volume definition to calculate the LF from galaxy 
datasets. The volume number density n obtained in equation ( fTSl l 
comes from the right hand side of Einstein's field equations and, 
therefore, it is written in terms of the proper volume. Thus, equa- 
tion ( fTsT i must have its volume units corrected by means of equa- 
tion ( |29l ) in order for the former to be correctly related to a se- 
lection function i// stemming from the LF. 

The relationship allowing us to calculate the completeness 
function used in this work may be written. 



ifr(z)^J(z)nc(z). 



(50) 



Now, if we go back to equation ( fT4b and the discussion in the 
paragraph below it, we can see that applying the completeness 
function by means of equation ( |49] | actually means replacing the 
theoretical comoving number density «c, based on a constant av- 
erage galaxy rest mass, with its observed counterpart ifr obtained 
directly from the LF. In that sense, all the observational quanti- 
ties obtained by applying J(z) to their theoretical counterparts 
will inherit that same empirical number count redshift evolu- 
tion encoded in the parametrization of the LF. As the goal of 
this paper is to consider the effect of redshift evolution on the 
relativistic number densities, such an approach should suffice. 
Also, in equation (l50b the completeness function is independent 



of volume units, therefore, if an observational quantity is ob- 
tained using J{z) by means of equation ( |49] l its original volume 
units dependence is preserved. 

The selection function may be written in terms of the LF, 



Hz) 



Jl(z) 



(f>{r)di, 



(51) 



where l(z) is the lower luminosity threshold below which the 
sources are not observed. The theoretical radial number density 
in comoving volume is given by. 



Vr 



3N(r) 



(52) 



Most observational quantities of interest studied in this paper 
require previous knowledge of the observed differential number 
counts [dN]„,„. Therefore, linking this quantity to its theoretical 
counterpart is an essential step in order to apply Ellis' equation 
view of this, to find [iiA^J^i.^ we identify T as dN in equa- 
tion (|49T l, yielding. 



[dNl„, = J{z)dN. 



(53) 



In essence this expression describes the same number counts as 
in equation ( ISOl l and both are in agreement with Ellis' (1971) key 
equation (flJl l. Considering now both equations dSOl l and (l53] l. it 
easily follows that. 



dN 
Iz 



- j{z)— - iML'^i 

dz ndz) dz 



(54) 



This is our key equation relating the relativistic theory to the 
observations. The number density of a proper volume is given 
by 



N_ 



(55) 



whereas the number density of a comoving volume is given by 
equation ( |52] |. Thus, we may rewrite equation (|54|) as follows. 



dN 
dz 



V^<p_dN_ 

Vp, n dz 



(56) 



It is important to point out that the volume transformation in this 
equation merely reflects the fact that J{z) is independent of vol- 
ume units, proper or comoving, and does not in itself change the 
volume units underlying both [dN/dz]„,„ and dN/dz when they 
are derived by means of Ellis' equation ( fT3l l. It is the number 
density n in equation ST3[ that defines the volume units of both 
[dN/dz]„i„ and dN/dz (see below). 

A further step towards writing these results in terms of 
the underlying cosmology is taken if we remember that RS03 
showed that Ellis' (1971) differential number counts ( fT3T l can be 
rewritten 



dN 



dy 



^ nid,y(l+z)dSl ^ . 
dz dz 

Thus, equation (l56l l becomes 



dN 




dz 


nhs 



-^{d,f{l+z)dQ. ^ 



Vr 



dz 



(57) 



(58) 



Note that in this expression [dN/dz]„i„ is given in terms of proper 
volume units because the number density n is written in terms of 
proper volume units. Also note that these expressions are gen- 
eral and independent of any cosmological model. The specific 
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cosmology will appear once we specialize the terms inside the 
brackets on the right hand side. 

Galaxy redshift surveys use bandwidth filters, and they 
sometimes include some sort of morphological classification. 
Therefore, a selection function tjj which considers a set of fil- 
ters W, and morphological types v, may be written as follows 
(see A07, eqs. 3-10), 



(59) 



Here f , is the fraction of galaxies in the dataset that were classi- 
fied with the morphological type v (see RS03, eq. 13), Al,, is the 
typical local rest-mass value for a galaxy of the morphological 
type V, and are constants introduced to avoid multiple count- 
ing of the same objects in the various filters, defined as follows 
(see RS03, eq. 46), 



Mz)=l, for 
and 

(z)<l, for W>1. 



(60) 



(61) 



Equation ( l6Tb simply states that when there exists more than one 
observed waveband (W > 1), gives the fraction of galax- 
ies in waveband W > 1 that are not counted in wavebands 
1,2,..., (W - 1). Considering these expressions, the complete- 
ness function defined in equation dSOl l can be rewritten 



«c(z) V 



(62) 



Substituting the expression (|59] | in equation (l58T l we obtain the 
following result. 



dN 




dz 





Vc 



dy 



-^id.rn + z)dn^ 



V, 



dz 



(63) 



This expression already appeared in A07. Here it has been 
reached through a more straightforward and simpler derivation 
and with the volume dependence explicitly shown. This means 
that volume conversions like the one discussed above in ^3.31 
may be needed if we are to calculate the observational differ- 
ential number count [dN/dz]„i,, in a consistent way. The above 
equation is shown here only to make explicit how the method- 
ology of previous work can be understood in relationship to that 
being developed in this paper In addition, it is important to note 
that equations ( l54l i. (l56b and ( l63T l are closely related. The first 
two are written in a more compact form, whereas the last one 
expands the morphological and bandwidth dependencies of the 
selection function and the relativistic features of the underlying 
four-dimensional spacetime, as indicated by the expression ( fSTl ). 

Equation ( l63T l also allows us to see how a theoretical mass 
evolution does not affect the differential number counts con- 
structed with LF parameters. This is based on the form of the 
summation over w, inasmuch as if we allow an evolution of the 
galactic mass by considering it dependent on the redshift, such 
a dependency will appear both in the numerator and the denom- 
inator of the summation and, therefore, cancels out, at least to 
first order However, although mass evolution does not directly 
alter equation (l63T l. it will indirectly change [dN/dz]„i,, as the LF 
parameters used in its construction will include some form of 
source evolution. See §2.2 of A07 for more details. 



4.2. Selection functions 

To calculate the completeness function and, as a consequence, 
other observational quantities, we need to compute the selection 
functions with respect to redshift in a given galaxy survey. G04 
and G06 fitted the LF Schechter parameters over the redshifts of 
5558 I-band selected galaxies in the FORS Deep Field dataset, 
photometrically measured down to an apparent magnitude limit 
of /,,B = 26.8. G04 and G06 also showed that the selection in the 
I-band is expected to miss less than 10% of the objects detected 
in the K-band, given that the AB-magnitudes of the I-band are 
half a magnitude deeper than those of the K-band, out to red- 
shift 6, beyond which no signal is detectable in the I-band due to 
the Lyman break. Heidt et al. (2001) also argue that the I-band 
selection minimizes biases like dust absorption. All galaxies in 
those studies were therefore selected in the I-band and then had 
their magnitudes for each of the five blue bands (1500 A, 2800 
A, m', g' and B) and the three red ones (r', /' and z') computed us- 
ing the best fitting SED given by their authors' photometric red- 
shift code convolved with the associated filter function. Gabash 
et al. (2004, 2006) determined the photometric redshifts by fit- 
ting template spectra to the measured fluxes on the optical and 
near infrared images of the galaxies. Heidt et al. (2001) reported 
that approximately 80% of the objects in that sample were clas- 
sified as Im/Irr-type, 15% as E/SO-type or Sa/Sb/Sc-type, and 
finally 5% as stars. 

The redshift ranges of 0.75 < z < 3.0 for the red bands and 
0.5 < z < 5.0 for the blue bands are large enough for checking 
possible observational inhomogeneities, since Rangel Lemos & 
Ribeiro (2008) showed that in the Einstein-de Sitter cosmology 
it is only theoretically possible, i.e., not observationally, to dis- 
tinguish these two features when z > 0.1, whereas for z > 0.5 
such a distinction becomes significant. Therefore, using galax- 
ies whose measured redshifts are mostly greater than z = 1 al- 
lows us to obtain data along the past light cone far enough from 
our present time hypersurface. One must remember that spa- 
tial homogeneity is defined on our constant time hypersurface 
whereas observational homogeneity occurs along our past light 
cone (Rangel Lemos & Ribeiro 2008; see also Ribeiro 1992, 
1995,2001,2005). 

To calculate the actual values of these functions, we begin 
by writing the result obtained in RS03 for the limited band- 
width version of the selection function of a given LF, fitted by a 
Schechter analytical profile, in terms of absolute magnitudes, 

rM'" (z) 

iA"(z) = 0.41nlO "" ^*(^)io°-41i+««««-«-*"'] X 

%J —CO 

xexp{-10°'^[*^'<'"'-'*'^VM"', (64) 

where, as discussed earlier, the index W indicates the bandwidth 
filter The equations for the redshift evolution of the LF parame- 
ters adopted by G04 and G06 are, 

fiz)^cf>;{'^+zf\ 

M*(z) = M*-H A'" ln(l +z), 
a(z) = a,„ 

with A* and B"' being the evolution parameters fitted for the dif- 
ferent W bands and M*, 0* and Uo the local (z « 0) values of 
the Schechter parameters as defined in G04 and G06. Since all 
galaxies were detected and selected in the I-band, we can write, 

Mliz) = M'Jx) = /,„„ - 5 log[d,(z)] -25+ A', (65) 

for a luminosity distance d^ given in Mpc. /,„„ is the limiting ap- 
parent magnitude of the I-band of the FDF survey and equals 



7 



A.S. Iribarrem et al.: Relativistic cosmology number densities and the LF 



to 26.8. Its reddening correction is A' — 0.035 (in Heidt et al. 
2001). We computed the selection functions for all eight bands of 
the dataset by means of simple numerical integrations at equally 
spaced values spanning the whole redshift interval. The errors 
were propagated quadratically, as mentioned before. The results 
are summarized in tables |2]to|4l Regarding the blue band dataset 
of G04, two distinct patterns of the selection functions in the dif- 
ferent bands are noticeable. The UV bands, 1500 A, 2800A, and 
u', evolve tightly with redshift, having values that are consistent 
with each other within the uncertainties, while at the same time 
assuming values outside the uncertainties of those in the blue 
optical bands g' and B. Therefore we chose to use the selection 
functions of the combined UV bands and those of the combined 
blue optical bands separately. The selection functions in the red- 
band dataset of G06, r', /' and z' are also combined. Once the 
combined selection functions have been obtained, we can calcu- 
late [dN/dz]„i,, by means of equations ( l54b and ( l62b . 

4.3. Observed number counts 

Following equation (l48l l. the best way to obtain an observed rel- 
ativistic number density for a given distance definition 
di is by calculating the observed number counts [N{z)]oi,s- These 
can be written as. 



[N(z)L 





dN' 


f 


dz' 



dz 



(66) 



The uncertainty of the cumulative number counts [A^]„,„ can be 
derived from the already determined uncertainty in the number 
counts [dN / dz]„i,,,. In this regard, the discussion about error anal- 
ysis in A07 (Appendix) showed that. 



where we can write. 



dN 




d^N 


-1 


dN 


dz 




dz^ 


6 

obs 


dz 



d^N 


d 


dN 


dz^ 


„„s dz 


dz 



d^N dJ dN 
dz} dz dz 



(67) 



(68) 



The derivative of the completeness function dJ{z)/dz can be ob- 
tained from equation (ISOl l and is proportional to the derivative of 
the selection function i/'"' itself with respect to redshift. This can 
be directly obtained from the parametrization of the LF. 

Next, we shall use the observed differential number counts 
[dN/dz]„i,, to study the relativistic radial distribution of the galax- 
ies in the G04 and G06 datasets in the same manner as in A07. 



5. Relativistic radial statistics of the FDF survey 

The observational values of the differential number counts 
[dN/dz]„i,, given in table|5]for all filters in the dataset permit us 
to evaluate the observational differential number densities [y]„i,,, 
for the various cosmological distance definitions in the FLRW 
model we are assuming. What it takes to do so is replacing 
dN/dz with [dN/dz]„i,,, in equation ( l44l i. which is exactly the 
same as applying the completeness function J(z) to the theo- 
retical values of the relativistic differential number densities y, . 
As discussed at the end of ^4.11 such application of the com- 
pleteness function naturally includes the observed evolution of 
the number counts in the selection functions of a given dataset. 

A look at equation (l44l i allows us to understand the sepa- 
rate roles of number counts and geometry, encoded in the term 



inside the brackets, in the relativistic differential densities y,. It 
shows that the underlying cosmological spacetime affects such 
quantities and that their theoretical departure from homogeneity 
at high redshifts is ultimately a consequence of the Universe's 
expansion, through the inclusion of the scale factor S in the dis- 
tance definitions, equations (fT9l l, ( l32l i. and (l34l i. The absence of 
the scale factor in the galaxy area distance da (eq. l33T l allows us 
to identify this distance with the radial coordinate r in this partic- 
ular spacetime. As a consequence, the differential number den- 
sities Jo obtained using this distance definition are not subject 
to the redshift evolution caused by expansion and will remain 
constant if a constant average galaxy rest mass is assumed. 

The fact that such differential densities depend on the ge- 
ometry of the Universe through the term inside the brackets 
of equation ( l44l i causes them to vary with the redshift differ- 
ently, and they do not necessarily remain constant. Such inho- 
mogeneities, however, are to be understood simply as a conse- 
quence of using relativistic distances, defined on the observer's 
past light cone. They are therefore easily reconciled with the 
Cosmological Principle, which requires large-scale spatial ho- 
mogeneity. As can be seen in figures 5 and 6 of A07, both differ- 
ential and integral number densities can be inhomogeneous, even 
when a constant average galaxy mass is assumed in an otherwise 
spatially homogeneous FLRW metric. Moreover, both quantities 
calculated using da remain homogeneous, as expected, since this 
distance in this particular spacetime is not affected by expansion, 
as discussed before. 

We can use this important result to further compare the rela- 
tivistic effects of considering densities down the past light cone 
in an expanding FLRW spacetime and the observational evolu- 
tion of the LF in that same hypersurface. The argument is simple: 
if we compute the ratio of the theoretical values of a given rel- 
ativistic density at a given redshift - say the differential number 
density using the luminosity distance d^ - to the one using dc, 
this indicates how much observational inhomogeneity is intro- 
duced through the use of d^ in the calculation of this relativistic 
density. Similarly, the ratio of the observational value of [ja],,!,, 
to its theoretical one, -yc,, indicates how much observational in- 
homogeneity the redshift evolution of the LF introduces. In fact, 
considering equations (l44l i and (|54] | we can see that such a ratio 
is simply the completeness function J(z) in a particular dataset, 
that is. 



[rc]„ 



[dN/dzlJAnjdc;)^ d{da)ldzY 
dN/dziMda^ d(da)/dz}-^ 
[dN/dzU, 



dN/dz 



= J(z). 



(69) 



Now, if there were no significant departure from homogeneity 
at a given redshift such ratios should be approximately unity. 
Therefore, if we subtract the values of those ratios at a given red- 
shift from one, this number should express how much inhomo- 
geneity that relativistic density shows at that redshift: the bigger 
the number, the bigger the inhomogeneity. Table |6] summarizes 
these results. At a fixed redshift entry one can compare the ef- 
fects on the homogeneity of the relativistic differential densities 
caused by the LF redshift evolution of the different filters in the 
G04 and G06 datasets. 

The results in Table |6] allow us to see that the evolution of 
the LF in the combined UV bands produces a bigger departure 
from homogeneity of the differential densities than the one in 
the combined optical bands, and this departure is even bigger 
than the one in the combined red ones. On the other hand, fixing 
a given ratio, say the observational-to-theoretical ratio for the 
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UV bands, [ycloi./Tc, one can observe that the departure from 
homogeneity of these relativistic differential densities is quite 
a common feature and tends to increase with redshift. Finally, 
by comparing the values of the purely theoretical jtlya ratio to 
the various observational-to-theoretical ones, one can investigate 
which contribution dominates at each redshift. That is particu- 
larly interesting, since the redshift evolution of any of the other 
differential densities can be obtained by a simple combination of 
those two separate effects, that of the expansion of the geome- 
try in the distance definition and that of the redshift evolution of 
the LF. For instance, let us consider the observational differential 
density in the combined UV filters calculated using the luminos- 
ity distance di, that is, [yiluv- Its departure from homogeneity 
can be investigated using the ratio [yt^uvlyG, as discussed above. 
Remembering that, by construction {ji^uv = Jm Jl, it follows, 

[yL]uv ^ Juvyi ^ [tgW Tt (70) 
Jg ya yo Tg' 

which is simply the combination of the observational-to- 
theoretical ratio for that dataset, namely the LF of the combined 
UV bands with the purely theoretical ratio for that distance def- 
inition, in that case, d^. The values in Table |6] indicate that the 
geometrical effect is comparable to the evolution of the LF in 
the whole redshift range of the G04 and G06 datasets. 

At this point, one should note that determining the spacetime 
geometry of the Universe or its density parameters is not the goal 
of the present paper We simply assume the ACDM metric, a 
necessary step in calculating the above mentioned geometrical 
terms, in order to see how the empirical evolution of the dif- 
ferential number counts [dN/dz]„i,, affects the theoretical results 
discussed above. 

Similarly, we can obtain the observational values of the inte- 
gral densities [y]*,,, by means of equation (l46b . or (l48b . because, 
by their very definitions, we have that [N]„i„ = Vc (/'(■z)- 

Figures[T]and|2]show graphs of the observational differential 
densities determined using the combined UV and optical bands 
in the G04 dataset plotted against the redshift, whereas figure 
[3] shows the results for those same densities versus the redshift 
in the combined red bands. The dependence of such relativis- 
tic densities on the distance definition, a known theoretical re- 
sult (Ribeiro 2005; A07), can be easily observed in these three 
figures. In addition, the dependence of [ydoi,, with the redshift 
in these same figures is solely due to the evolution of the mass 
function, since its relativistic volume definition remains constant 
in the assumed FLRW spacetime. We note that values of [yo]„i,., 
obtained from the LF in the combined red bands (figure|3]l seem 
less dependent on the redshift than those in the UV and optical 
bands (figures [1] and |2]i. 

Similar graphs for the observational integral densities in the 
UV, optical and red combined bands versus the redshift are re- 
spectively shown in figures |4l|5] and|6] where their dependency 
with the distance definition can also be clearly seen. Although 
less pronounced in these smoothed out relativistic densities, the 
effect of the evolution of the mass function can still be observed 
in these three graphs, as [y*j]„i,,, is also unaffected by the geomet- 
rical effects of considering the densities on the past null light 
cone of an FLRW spacetime. In other words, the LF redshift 
evolution seems to further enhance the inhomogeneity of these 
densities, making observationally inhomogeneous even the the- 
oretically homogeneous yc and y*. We observe that in figure |6l 
[TgIo'" is also less dependent on the redshift than in the UV and 
optical bands shown in figures |4] and |5] 

Figures |7] [8] and|9]plot the observational differential and in- 
tegral densities built with the luminosity distance in the com- 



bined UV, optical and red bands, respectively, versus the lumi- 
nosity distance. A power-law pattern for both densities in the 
three combined bands can be observed when one compares the 
data points with the solid line drawn just as reference in all three 
plots. The slope of the reference lines is the same in all bands, 
which indicates that there is little change in the slope of the data 
points. This possibly indicates that such an effect is predomi- 
nantly relativistic. 

A similar pattern can be found in figures [TOl [TT] and [12] 
which plot the observational differential and integral densities 
using the redshift distance against that same distance in the com- 
bined UV, optical and red bands respectively. We found a similar 
indication of a power-law pattern in both densities at high red- 
shifts, but the actual power exponent of the distributions is differ- 
ent from the one using the luminosity distance. This comes from 
the fact that these two relativistic distances are affected differ- 
ently by the expansion of the FLRW spacetime. Comparing the 
results shown in figure [TT] with the ones in the UV band in fig- 
ure[TO]we can see that the slopes of the data points should agree 
with one another. Similarly as in figures [10] and [TT] the power 
law exponent shown in figure [12] seems to be the same for both 
[Tz],*i and [y^],,,,, in the given distance definition. 

One should notice that the slope of the straight lines used 
as references for the data points in figures [7] to [12] do not vary 
drastically over all bands of the FDF datasets, both for luminos- 
ity and redshift distance definitions. This possibly indicates that 
such an effect is predominantly relativistic. 

The differential density constructed with the area distance t/^ 
becomes discontinuous atz~ 1.5 due to the fact that, by defini- 
tion, d,i has a maximum at that redshift and, therefore, its deriva- 
tive with respect to z becomes zero, yielding an undefined [y^l^tj 
at high redshifts. In addition, differently from the decrease of y* 
and 7* at higher redshifts, the integral density y* increases with 
z. Such pathologies have already been previously detected (see 
figs. 5 and 6 of A07). They seem to render the number densities 
defined by d,i as unphysical. 

The differential and integral densities defined by dc lack the 
geometrical effect of expansion in the particular spacetime un- 
der consideration in this paper and, therefore, are understood to 
be unsuitable if one wants to fully characterize the relativistic 
power-law patterns arising from the combination of both space- 
time expansion and LF evolution. 

Finally, it is apparent that the differential densities [yi]„,„ and 
[y-]„,„ do not appear to exhibit a power-law decay as the inte- 
gral ones do, especially at redshift values nearing the sample's 
depth. This is very similar to what found in A07. This behavior 
is possibly a result of noisier data at the redshift limit of the sam- 
ples. By definition, the differential densities measure the rate of 
growth in number counts, rendering them more sensitive to local 
fluctuations, whereas the integral densities indicate the change in 
number counts for entire observational volumes, rendering them 
less sensitive to the same fluctuations. 

The indication of power-law behavior may be the result of 
observational, not necessarily spatial, inhomogeneity as we look 
down our past light cone, combined with the incompleteness of 
galaxy counts at higher redshifts. It may also be a result of other 
causes which we are investigating. 

6. Conclusion 

In this paper we have detailed a generalization of the framework 
connecting the relativistic cosmology number count theory with 
the astronomical data extracted from the galaxy luminosity func- 
tion (LF), as proposed by Ribeiro & Stoeger (2003). This frame- 
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work was used by Albani et al. (2007) to extract the observed 
differential number counts from the LF data and to study the 
observational inhomogeneities in the relativistic radial densities 
using different distance definitions. Here we have carried out an 
analysis of number counts using the LF dataset of Gabasch et al. 
(2004) in the redshift range 0.5 < z < 5.0, and of Gabasch et 
al. (2006) in the redshift range 0.75 < z < 3.0. In doing so we 
focused on the observational inhomogeneities in the relativistic 
radial statistics of the distribution of the galaxies by means of 
the differential densities y, (/ = a, c, l, z) and integral densi- 
ties y", both using the various cosmological distance definitions 
di, as well as the selection functions and the differential number 
counts [dN/dz]„bs obtained from the LF of G04 and G06. 

We confirmed the dependence of such empirically based rel- 
ativistic densities on the distance definition, a geometrical con- 
sequence of the expansion of the Universe, which also is the 
cause of the observational inhomogeneities present even in the- 
oretical analysis of the spatially homogeneous FLRW spacetime. 
Furthermore, as expected, the effect of the redshift evolution of 
the LF seems to increase the observational inhomogeneity of 
these densities. 

We found evidence of a power-law pattern in the behavior 
of [yi]„i, and [y;]„,„, relative to their respective cosmological dis- 
tance measures d^ and da, similar to that previously detected in 
A07. However, this pattern is somewhat less pronounced than 
the results of A07. This difference could be due to the use of a 
different, and somewhat more complete, redshift samples used 
in this paper, especially the red galaxies of G06. 

Finally, we should note the dependence of these results 
and their interpretation on the assumed FLRW cosmology. 
Therefore, we look forward to finding out what an analysis 
of this data with reference to an inhomogeneous cosmological 
model (e.g., a Lemaitre-Tolman-Bondi model) would reveal. If 
the universe is not almost-FLRW or if the data are not probing 
almost-FLRW scales, then hidden in the results are the effects 
of having chosen the wrong cosmological model. In addition, 
what we interpret as LF evolution, which is the same as number 
evolution, but with some indirect luminosity evolution effects, 
undoubtedly also reflects a lot of unobserved, and therefore un- 
counted, galaxies. Strictly speaking, if the differential and in- 
tegral densities are comoving densities, then jo and y* should 
increase with redshift, as there should be more galaxies in a co- 
moving volume in the past that there are now due to mergers 
(more, but smaller, galaxies in the past). So we are obviously 
not observing as many galaxies as we go farther out. Of course, 
that would mean that the average mass per galaxy decreases with 
redshift, an effect that should be implicit in the empirical LF, 
but which requires further theoretical considerations to be made 
explicit in our model. Including these points in our general ap- 
proach to the problem of connecting relativistic cosmology the- 
ory and the empirically determined LF is indeed necessary in or- 
der to make the present study even more realistic. We intend to 
deal with these issues, as well as others like the possible meaning 
of power law patterns in the differential and integral densities, in 
forthcoming papers. 
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Appendix A: Numerical algorithm 

This appendix describes step by step how to apply the theoretical 
and numerical methodologies described in sections ^ ^and 
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Table 2. Selection functions for the FDF blue UV bands 



Table 4. Selection functions for the FDF red bands 
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Table 3. Selection functions for FDF blue optical bands 



redshift 
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^in order to obtain the results used in this paper. Steps 1 to 6 
involve the calculation of the theoretical quantities, as discussed 
in ^and ^ Steps 7 to 10 show how to obtain the corresponding 
observational quantities for a given LF, as discussed in ^ 

Step 1 - Coupled differential equations. Obtain the scale fac- 
tor S (r) and the cumulative number count N(r) as a function of 
the radial coordinate, or comoving distance r, by solving the cou- 
pled equations ( l20l l and (|23] |. We used the fourth order Runge- 
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Table 5. Differential number counts 
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Table 6. Spacetime expansion vs. galaxy evolution effects on 
relativistic homogeneity 
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Kutta method to obtain values of S (r) and N(r) at fixed r incre- 
ments, assuming as initial conditions A^(0) - and 5 (0) = 1 ; 

Step 2 -Auxiliary quantities. Obtain the redshift z(r) for each 
increment of r using equation (l24l) and the corresponding value 
of the scale factor S(r). Do the same for the area distance d^ir) 
using equation ST% . the differential number count dN/dz(r) us- 
ing equation ( |28] ), and the comoving volume number density 
ndr) using equation ( l52] l: 

Step 3 - Cosmological distances. Having the values of z(r) 
and d^ir) calculated for each r increment, obtain the values of 
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Fig. 1. Observational relativistic differential densities versus red- 
shift in the combined UV bands of the FDF dataset of G04. 
Symbols are as in the legend. 



Fig. 3. Observational relativistic differential densities versus red- 
shift in the combined red bands of the FDF dataset of G06. 
Symbols are as in the legend. 
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Fig. 2. Observational relativistic differential densities versus red- 
shift in the combined optical bands of the FDF dataset of G04. 
Symbols are as in the legend. 



Fig. 4. Observational relativistic integral densities versus red- 
shift in the combined UV bands of the FDF dataset of G04. 
Symbols are as in the legend. 



the galaxy area distance do{r), and the luminosity distance dtir) 
by direct application of equation ( l30t . Etherington's reciprocity 
law itself. The values of the redshift distance d,(r) can be easily 
obtained for each r increment using equation ( |3TI ): 

Step 4 - Derivatives of the distances with respect to redshift. 
For each r increment, calculate the values of the derivatives with 
respect to redshift of the distances obtained in the previous step 
by means of equations ( |36] |. ( [39] l and ( l40l l: 

Step 5 - Theoretical differential densities. Obtain the values 
of the theoretical differential densities y,{r) for each of the dis- 
tance definitions d, (; = a, g, l, z) calculated previously for each 
r increment, by means of equation (l44l l: 

Step 6 - Theoretical integral densities. Calculate for each 
distance d^{r) entry its corresponding spherical volume V^r), by 
means of equations (l43b . then use those values, together with 



the ones for the cumulative number count N{r), to obtain the 
values of the integrated differential densities y*{r), by means of 
equation ( l48T i: 

Step 7 - Selection functions. Select a bin containing a group 
of redshift values in a certain interval {Za^b\, from a point Za to 
a point Zb, taken from the ones obtained in step 2, and obtain 
the corresponding values of the selection functions i//^iza^h) for 
each filter W using equation ( |64] |. with the appropriate values for 
the LF's Schechter parameters (f>*(Za^b), M*{za^b), and a(za^h), 
and different absolute magnitude limits MJ^^Jza^t), given by 
equation ( l65l i: 

Step 8 - Completeness functions. Calculate the values of the 
completeness function J^iza^b) by means of equation (ISOl l us- 
ing the selection functions tfr'^iza^t) calculated in step 7 and the 
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Fig. 5. Observational relativistic integral densities versus red- 
shift in the combined optical bands of the FDF dataset of G04. 
Symbols are as in the legend. 



10^ 



10" 



- 10-^ 

I 

u 
Q. 

SlO^ 



^10-^ 



10' 



10 



i i 



° ° [7.4 1 red 

• • baled 

' ' [jz]red 



^ i i i i 

' f 1 . , 



0.7 



1.0 



2.0 



3.0 



Fig. 6. Observational relativistic integral densities versus red- 
shift in the combined red bands of the FDF dataset of G06. 
Symbols are as in the legend. 



corresponding values of the comoving volume number density 
nc[Za^bif)] obtained in step 2; 

Step 9 - Differential and cumulative number counts. Obtain 
the observational differential number count [dN/dz]^^X^a^b) in 
each filter W by means of equation ( |54] | and the observational cu- 
mulative number count [A^l^Xzo^/,), by means of equation ( l66] l: 

Step 10 - Observational differential and integral densities. 
Obtain the observational differential densities {y^'^,^,(Za^h), and 
the integrated differential densities [y'lJ^X^a^*)' using the 
values of [dN I dzYJza^h) and [A^],l(Za^i) in place of their the- 
oretical counterparts dN/dz and in equations (l44l i and ( |48]) . 
respectively. 
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Fig. 7. Observational relativistic differential and integral densi- 
ties versus luminosity distance in the combined UV bands of the 
FDF dataset of G04. The solid line was drawn just as reference, 
to show the apparent power-law behavior of both relativistic den- 
sities in that distance definition at high redshifts. Symbols are as 
in the legend. 




Fig. 8. Observational relativistic differential and integral densi- 
ties versus luminosity distance in the combined optical bands of 
the FDF dataset of G04. The solid line is plotted just as a refer- 
ence to the power law behavior stated in the legend. Symbols are 
as in the legend. 
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Fig. 9. Observational relativistic densities versus luminosity dis- 
tance in the combined red bands of the FDF dataset of G06. 
The solid line is shown as a reference to the power law behavior 
stated in the legend. Symbols are as in the legend. 
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Fig. 11. Observational relativistic differential and integral densi- 
ties versus redshift distance in the combined optical bands of the 
FDF dataset of G04. The solid line is shown as a reference to the 
power law behavior stated in the legend. Symbols are as in the 
legend. 
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Fig. 10. Observational relativistic differential and integral den- 
sities versus redshift distance in the combined UV bands of the 
FDF dataset of G04. The solid Une is shown just as a reference 
to the power law behavior stated in the legend. Symbols are as 
in the legend. 
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Fig. 12. Observational relativistic differential and integral den- 
sities versus redshift distance in the combined red bands of the 
FDF dataset of G06. The solid line is plotted as a reference to 
the power law behavior stated in the legend. Symbols are as in 
the legend. 
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